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In quantum adiabatic algorithm, as the adiabatic parameter s(t ) changes slowly from zero to one 
with finite rate, a transition to excited states inevitably occurs and this induces an intrinsic compu¬ 
tational error. We show that this computational error depends not only on the total computation 
time T but also on the time derivatives of the adiabatic parameter s(t) at the beginning and the 
end of evolution. Previous work (Phys. Rev. A 82, 052305) also suggested this result. With six 
typical paths, we systematically demonstrate how to optimally design an adiabatic path to reduce 
the computational errors. Our method has a clear physical picture and also explains the pattern of 
computational error. In this paper we focus on quantum adiabatic search algorithm although our 
results are general. 


PACS numbers: 03.67.Ac, 03.67.Lx, 89.70.Eg 

I. INTRODUCTION 

Quantum adiabatic algorithm was proposed in 2000 
by Farhi et aZ.jTj as an alternative paradigm of quan¬ 
tum computing to quantum circuit algorithm [2]. It 
works by constructing a time-dependent Hamiltonian 
that evolves slowly from the initial Hamiltonian to the 
problem Hamiltonian. The ground state of the initial 
Hamiltonian is easy to find and the answer of the intent 
problem is encoded in the problem Hamiltonian. It is 
ensured by the quantum adiabatic theorem that if the 
Hamiltonian evolves slowly enough, the system will stay 
in the ground state and evolve into the ground state of 
the problem Hamiltonian. When it is applied to the 
search problem, the quantum adiabatic algorithm has 
been shown to be 0(y/~N) [31 0 ], which is as powerful 
as Grover’s algorithm [5; and is quadratic speedup over 
the classical search algorithm. In general, the adiabatic 
algorithm has been shown to have the potential to solve 
NP-hard problem El Cl¬ 
in addition to its speed up against classical comput¬ 
ing, the quantum adiabatic algorithm also has capacity 
to remain robust against environment noise [8]. Some 
practical architectures for quantum adiabatic algorithms 
were proposed [SJ; in 2013 D-Wave company claimed that 
they built a quantum computer based on quantum adia¬ 
batic algorithm ma¬ 
in quantum adiabatic algorithm, there are two types 
of computational errors. One is the extrinsic error, which 
is caused by the environment. The other is the intrinsic 
error: as the algorithm has to be run in a finite compu¬ 
tation time T, the adiabatic Hamiltonian must change 
in a finite rate and this inevitably will induce transition 
to excited states and cause computational error. The in¬ 
trinsic computational error depends entirely on how the 
adiabatic path s(t) is chosen. The most popular choice so 
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far is the linear path, s(t ) = t/T. Other choices were pro¬ 
posed in literature m People has also tried to optimize 
the adiabatic path s(t) using geometrization [4]. 

According to the hierarchical theory of quantum adia¬ 
batic evolution [12], the intrinsic error depends crucially 
on the time derivatives of s(t) at the beginning and the 
end of the evolution. This fact was also pointed out in 
Ref- [H. E3- In this work we choose six typical adi¬ 
abatic paths s(t) to systematically demonstrate how to 
optimally control these time derivatives to reduce the 
computational error. Also we give a physical picture 
about the origin of computational error’s oscillation. We 
focus on the quantum adiabatic search algorithm and nu¬ 
merically compare the computational errors for different 
adiabatic paths. We find that the cubic path is the best 
among the six chosen path and it can reduce the error 
by orders of magnitude compared to the popular linear 
path and sinusoidal path. 

Our paper is organized as follows: for the sake of self¬ 
containment, we first briefly introduce the quantum adi¬ 
abatic algorithm, in particular, the quantum adiabatic 
search algorithm, and the hierarchical theory of quan¬ 
tum adiabatic evolution, respectively, in Sections II and 
III. In Section IV, we show that the Hamiltonian for the 
quantum adiabatic search algorithm can be reduced to a 
spin-1/2 Hamiltonian and apply the hierarchical theory. 
In Section V, six typical adiabatic paths are chosen and 
they are categorized into three groups. Our main numer¬ 
ical results are shown in Section VI. We finally conclude 
in Section VII. 


II. QUANTUM ADIABATIC SEARCH 
ALGORITHM 

The adiabatic quantum computation was first intro¬ 
duced in 2000 [1] based on the quantum adiabatic theo¬ 
rem |14| . In a quantum adiabatic algorithm, the solution 
of a problem is encoded into the ground state of the prob¬ 
lem Hamiltonian H p . An initial Hamiltonian H is chosen 
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so that its ground state can be easily found and set up. 
The total adiabatic Hamiltonian is constructed by link¬ 
ing the initial Hamiltonian with the problem Hamiltonian 
with a path as follows 


H s (t) = (l-s(t))H b + s(t)H p , (1) 

where s(t) £ [0,1] is the adiabatic parameter. The sys¬ 
tem is prepared in the ground state of Hb. As s changes 
slowly from zero to one, the quantum adiabatic theorem 
ensures that the system stays in the ground state of H s 
and eventually arrives at the ground state of H p , the so¬ 
lution. 

In the search problem, the task is to locate M 
marked items out of N randomly arranged items. On 
a quantum computer, we use a set of orthonormal ba¬ 
sis 11), 12),..., \N) to denote the N unsorted items. The 
problem Hamiltonian can be constructed as [3] 

H P = I” Y l m )H; ( 2 ) 

mGA1 


where M. is the set of marked items. This Hamiltonian 
is the projection operator to a subspace orthogonal to 
the subspace spanned by its ground state 

can be any state from subspace If we choose 

the initial state to be equally contributed from the or¬ 


thonormal basis as 


\ipo) = Eili |i), then the initial 


Hamiltonian is 


H b = 1 - \ip 0 )(Po\ = 1 - (3) 

hj 


Thus the quantum adiabatic Hamiltonian for search is 

h s = i - 1 ^ Y l»)01 - 5 (t) Y l m )H • ( 4 ) 

i,j m£A4 


Note that the unit of Hamiltonian which rules the 
quantum computation depends on what kind of system 
is used for realization. Suppose a system has a charac¬ 
teristic time r; the unit of Hamiltonian is H/t. For con¬ 
venience, in the following derivation, we set r = K = 1. 

In the above Hamiltonian, if s(t) changes from zero to 
one infinitely slowly, the system will stay strictly in the 
ground state and the solution can be found without any 
error if there is no environment noise. This is dictated 
by the quantum adiabatic theorem. However, we want to 
know the solution as fast as possible. This means that s 
has to change from zero to one in a finite time T, causing 
a small transition to excited states. At the end, there is 
an inevitable error in the solution. For the search prob¬ 
lem, we define the intrinsic computational error as m 

5=l~(m\P\ 1 p(T)), (5) 


where the projection operator P is defined as 


P= Y 

m^M. 


( 6 ) 


The main purpose of this work is to reduce the com¬ 
putational error S by optimizing the adiabatic path s(t). 
The simplest and also the most popular choice is the lin¬ 
ear path s = t/T. This is almost the worst among the 
easy choices, which includes sinusoidal path. According 
to the newly developed hierarchical theory of quantum 
adiabatic evolution |12] , the error 6 depends crucially 
how the time derivatives of s(t) at the beginning and 
the end of the adiabatic evolution. In this work, several 
adiabatic paths s(t) are designed and the errors caused 
by these paths are computed and compared to the error 
by the linear path. The error can be reduced by orders 
of magnitude. 


III. HIERARCHICAL THEORY OF QUANTUM 
ADIABATIC EVOLUTION 



FIG. 1: (color online)Adiabatic evolution trajectories in the 
projective Hilbert space. The red lines are the zeroth-order 
trajectories, that is, the trajectories follow strictly the quan¬ 
tum adiabatic theorem. The blue lines are the first-order 
trajectories, which are shifted from the zeroth-order ones by 
a small amount proportional to s. The black lines are the 
second-order trajectories. When s is not zero at t = 0, the sys¬ 
tem will oscillate around the first-order trajectory (the green 
line in (a)). Similarly, when s is not zero at t = 0, the system 
will oscillate around the second-order trajectory (the green 
line in (b)). 


The quantum adiabatic theorem was proved in 1928 by 
Born and Fock m- This theorem ensures that a system 
starting in the ground state will stay in the ground state 
when the adiabatic parameter s changes slowly. How¬ 
ever, this is mathematically true only when the changing 
rate of s is infinitesimally small. In any practical situ- 
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ation, for example, quantum adiabatic computing, the 
adiabatic parameter s has to change with a small but 
finite rate, this will cause a small transition to excited 
states, resulting a deviation from the quantum adiabatic 
theorem. In Ref. [T2], a hierarchical theory was devel¬ 
oped to compute the deviation order by order. 

The results are schematically illustrated in Fig. |T]with 
trajectories in the projective Hilbert space (the overall 
phase is not important) [Tbfl 17[ . At the zeroth order, the 
system follows the trajectory dictated by the quantum 
adiabatic theorem. At the first order, the system oscil¬ 
lates with a small amplitude around a trajectory slightly 
shifted from the zeroth-order trajectory. The small shift 
is proportional to s, the first-order time derivative of 
s(t); the small oscillating amplitude is determined by s 
at t = 0. This is illustrated in Fig. 3 a )- 

To reduce and also to better control the computational 
error S, we can set s = 0 at the beginning. In this 
case, there are no oscillations and the system will fol¬ 
low a smooth first-order trajectory. As the shift between 
the first-order trajectory and the zeroth-order trajectory 
is proportional to s, we can reduce the error S in the 
first order to zero by choosing a s(t) such that s = 0 at 
t = T. For such a path s(t), the error S is of the second- 
order, determined by s, the second-order time derivative 
of s(t). To further reduce the error, we can repeat the 
above procedure by choosing a path s(t) such that s = 0 
at t = 0, T. This strategy is very successful as we shall 
see in the following sections. 

We note here that the crucial role played by the time 
derivatives was also pointed out in Ref. mm3 with a 
different approach. 


where r = M/N. In other words, the search Hamiltonian 
H s is reduced to a spin-1/2 Hamiltonian H s . Fig. [ 2 ] 
shows the eigenvalues of both Hamiltonians, H s and H s . 
It is clear that two eigenvalues of H s (diamonds) are 
identical to two of the eigenvalues of H s (red solid lines). 
Due to the permutation symmetry, all other eigen¬ 
states of H s (non-red solid lines) do not participate in 
the dynamical evolution when s changes from zero to one. 



FIG. 2: (color online)Eigenvalues of the search Hamiltonian, 
(a) N = 100, M = 1; (b) N = 100, M — 2. The solid lines are 
for H a and the diamonds for H s . The time unit depends on 
the system used to realize our algorithm. Suppose the system 
has a characteristic time r, then the unit of energy is h/r. In 
this article, we have set h = r = 1. 


IV. REDUCED SEARCH HAMILTONIAN AND 
ANALYTICAL RESULTS 


Let us come back to the search Hamiltonian and see 
how the hierarchical theory outlined in the last section 
can be applied successfully to this problem. Because of 
the permutation symmetry of the search Hamiltonian H s , 
for the given initial state, the quantum state at any given 
time has the form [18] 




E i u > + -fe E m ■ ( 7 ) 


y/N-M 


u£M. 


Vm 


m£Ai 


To apply the hierarchical theory in Ref. [12], we refor¬ 
mulate our problem in the projective Hilbert space [13- 
H7] . We rewrite the state as 


1 ^ 2 ) = 




( 10 ) 


and define p = arg(A 2 ) — arg(Ai) £ [0,27 t], q = |A 2 | 2 £ 
[0,1]. The dynamics of our system in the projective 
Hilbert space can be completely determined by p and q. 
In terms of p and g, the classical Hamiltonian becomes 


The Schrodinger equation governing the search algorithm 
becomes 




dt u, 



( 8 ) 


H c = ^ + y/q( 1 — q) cos(p)(2r + 2s — 2 rs — 1) 

-2y/q(l - q) sin(p) y/r( 1 - r)(l - s) , (11) 


with 

ft = ( r(l - s(t )) + s s(t))\ 

\-y/r(l - r)(l - s(t)) (1 — r )(l — s (t)) ) 

(9) 


where H c = (^ 2 |iJ s |^ 2 ). 


The ground state of the reduced search Hamiltonian 
is a fixed point in the projective Hilbert space where the 
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overall phase is removed. The fixed point is given by 


2. Sinusoidal path 


<1 = o 


7T — arctan( 


2>/(! - r)r(l - s) 


2r + 2s — 2rs — 1 


)■ (s> 


1 - 2r\ 


2 — 2r 


P = 


arctan( 


2y / ( 1 ~ r)r(l - s) 
2 rs — 2r — 2s + 1 


)■ (*< 


1 - 2r 


2-2 r 

( 12 ) 

where we have assumed that r <C 1. When s changes 
infinitesimally slowly, this fixed point traverses the adia¬ 
batic trajectory (the red lines in Fig. l]). When s changes 
with a finite but small rate, the actual dynamics will devi¬ 
ate from the adiabatic trajectory of the fixed point. The 
averaged first-order and second-order deviations are j!2j 


1. First Order Deviation 


[ 0 t < 0 

s 2 (t) = ^ sin(f£) 0 < t < T (17) 

1 t>T 


3. Square path 


s 3 {t) 


0 t < 0 

3(A) 2 -2(A) 3 0 <t<T 

1 t>T 


(18) 


4. Sinusoidal square path 


[0 t < 0 

s i(t) = < sin 2 ^) 0<t<T 
I 1 t>T 


(19) 



2. Second Order Deviation 


(13) 


5. Sinusoidal cubic path 


[0 t < 0 

s 5 (t) = < sin 3 (^p) 0 <t<T 
I 1 t>T 


( 20 ) 


^ = i2 x MT^) 3 ’ ,r ~ 1)( y 4,) ~ 

where A = 1 + 4(r — l)s — 4(r — l)s 2 . Ai is the ith 
order deviation of p, and Bi is the ith order deviation 
of q. Detailed derivation of ( Ai,Bi ) can be found in the 
appendix. The first-order deviation as depicted by the 
blue lines in Fig. |Tj is proportional to s. The second- 
order deviation is depicted by the black line in Fig. [ljb) 
and it is proportional to s if s is zero. 



V. ADIABATIC PATHS 


6. Cubic path 

ro t< o 

s 6 (t) = < 6(|0 5 - 15(^) 4 + 10(T) 3 0 <t<T (21) 

[1 t>T 

Among these six paths, linear path and sinusoidal path 
are of zeroth order; the first-order paths are square path, 
sinusoidal square path, and sinusoidal cubic path; cubic 
path is second order. As we shall see, the higher-order 
path can lead to smaller computational error by orders 
of magnitude. Note that only linear path and sinusoidal 
path of those six paths were studied before. mi 


We will show in the last two sections that the time 
derivatives of s at the beginning and end play a crucial 
role in determining the computational errors. We define 
these time derivatives as 

d n <s d n s 

c " = ^ (0 > ■ < 15 > 

With c n and d n , we categorize the evolution paths in the 
following ways: the path is called nth-order if all of its c m 
and d m for m < n are zero while either c n or d n are not 
zero. The zeroth-order path has either c\ 0 or d\ ^ 0. 
We propose the following six adiabatic paths to illustrate 
our ideas: 

1. Linear path 

f 0 t < 0 

ai(*)=< | 0 <t<T (16) 

1 t > T 


VI. NUMERICAL RESULTS 

The Schrodinger equation ([8]) is solved numerically for 
the six adiabatic paths for various values of T. The re¬ 
sults are analyzed, explained, and compared in this sec¬ 
tion. They convincingly show that higher-order paths 
can reduce computational errors by orders of magnitude. 

Fig- 0 shows how the deviation from the ground state 
(fixed point) changes with time for linear path and si¬ 
nusoidal path. It is clear from this figure that the devi¬ 
ations oscillate around the first-order analytical results. 
And the amplitude of this oscillation is conserved, which 
is an result of conservation of action[T9T. This is expected 
from the hierarchical theory (see Fig. [l])P3 as both lin¬ 
ear path and sinusoidal path are zeroth-order path with 
nonzero C\. It is interesting to compare the results for 
linear path and sinusoidal path. The deviation for the si¬ 
nusoidal path is larger than that for the linear path in the 
middle of the evolution; however, the computational error 
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FIG. 3: (color online) Deviations from the ground state as a 
function of time for the linear path and the sinusoidal path. 
The analytical results for the first-order deviation are plotted 
for comparison. (N = 10, M = 1, T = 1000) 


(i.e., the deviation at t = T) appears slightly smaller for 
sinusoidal path. The reason is that s decreases smoothly 
toward zero for the sinusoidal path. Although the reduc¬ 
tion of the computational error is not much, it already 
shows the possibility to reduce the computational error 
by optimally designing s(t) for a given T. This kind of 
reduction can be of order of magnitude when we choose 
a path of higher-order as we shall see next. 

We have designed three first-order paths, square path, 
sinusoidal square path, and sinusoidal cubic path whose 
ci and d\ are zero. As an example, our numerical results 
for square path is plotted in Figj4j The results in Fig0a) 
are deviation from the ground state. As the first-order 
derivative ci and d\ are zero, the deviation is of second 
order and it oscillates around the analytical second-order 
deviation (see the inset of Figjija)). To see these oscilla¬ 
tions more clearly, we have plotted the difference between 
our numerical results and the analytical second-order de¬ 
viation in Figjdjb)), where the oscillation pattern is seen 
to have a kink in the middle of the evolution. Since 
the second order derivatives C 2 and are very small, 
the deviation is very small, order of magnitude smaller 
than the ones in Fig. [3j The deviation pattern for sinu¬ 
soidal square path and sinusoidal cubic path are similar 
to square path. To avoid the overcrowding of figures, we 
only plotted the results for square path. 

We finally look at the results for cubic path, which is 
the only second-order path among the six. From the dis¬ 
cussion above, we know that the main deviation of cubic 
path is of the third order because both its C 2 and are 
zero. Therefore, the deviation of cubic path should be 
very different from square path. For comparison, the nu- 



FIG. 4: (color online) (a) Deviation from the ground state as 
a function of time for the square path. The analytical results 
for the second-order deviation are plotted for comparison, (b) 
Difference between numerical result of p and its second order 
analytical result as a function of time. (N = 10, M = 1, 
T = 1000) 


merical results for cubic path are plotted in Figj5] in a 
similar fashion as for square path. There are still oscil¬ 
lations around the second-order analytical result as seen 
from the inset of Fig(5ja) . However, the oscillations have 
a very different pattern: as C 2 is zero, the oscillation 
starts at zero with a much smaller amplitude and fre¬ 
quency. In Fig(5jb), where the difference is plotted, the 
oscillations with small oscillations are not even visible 
due to a large kink around t = 500. This is very differ¬ 
ent from FigQb)). Overall, we see that the deviation of 
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FIG. 5: (color online) (a) Deviation from the ground state as 
a function of time for the square path. The analytical results 
for the second-order deviation are plotted for comparison, (b) 
Difference between numerical result of p and its second order 
analytical result as a function of time. (N = 10, M = 1, 
T = 1000) 


FIG. 6: (color online) The computation error as a function 
of computation time T for the linear path and the sinusoidal 
path. (N = 100, M = 1) 



cubic path is much smaller than square path. 

For a given total computation time T, the success of 
an algorithm depends on the computational error. The 
smaller the error the better the algorithm. The compu¬ 
tational errors for all the six proposed paths are com¬ 
puted and plotted against the computation time T in 
Figs. [B0&|5| Fig. U shows the results for the two path 
of the zeroth order. The results for the three first-order 
paths are shown in Fig [7] The cubic path is of the second 
order and has the smallest computational error as seen 
in Fig. [8j The oscillations in these three figure are orig¬ 
inated from the oscillations in Figs. [3]&|4j they are not 


FIG. 7: (color online) The computation error as a function of 
computation time T for the square path, the sinusoidal square 
path, and the sinusoidal cubic path. (N = 100, M = 1) 


essential. 

We can view the results in Figs. |6|7[ &|8| in a different 
angle: for a given allowed computational error, which 
path has the shortest computation time? For this angle, 
we have averaged out these oscillations and combined 
the results in Figs. |6|7[ tf8l into Fig. [9] In Fig. |9j we see 
clearly that the relation between the computation time 
T and the computational errors differs slightly for paths 
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FIG. 8: (color online) The computation error as a function of 
computation time T for the cubic path. (N = 100, M = 1) 



FIG. 9: (color online)The smoothed computational error with 
total evolution time of six paths. (N=100, M=l) 

of the same order. However, for path of different orders, 
the relation is very different: for a certain allowed com¬ 
putational error, the computation time T can differ by 
orders of magnitude; similarly, for a given computation 


time T, the error can differ by orders of magnitude. For 
example, at T = 100, the computational error for cubic 
path is almost seven orders of magnitude smaller than 
the popular linear path. If we want to limit the com¬ 
putational error to 10 -6 , it only takes T = 25 for cubic 
path; for sinusoidal square path, T is around 46; for the 
popular linear path, it would take a much longer time 
T = 650. 

Although the oscillations are not essential, they do ex¬ 
hibit interesting patterns. The computational errors for 
four paths, the linear path, the square path, the sinu¬ 
soidal square path, and the cubic path, have very similar 
oscillation patterns while the other two paths, the sinu¬ 
soidal path and the sinusoidal cubic path, share another 
oscillation pattern. The pattern is largely depend on the 
relationship between \c n \ and |d n |, the derivatives of s 
at the beginning and end. Take zeroth-order path as an 
example. If |ci| = |di|, the actual trajectory can inter¬ 
sect the adiabatic trajectory for some values of T. At 
those T, the computational error is extremely small, as 
shown in Fig. 1(a). However, if |ci| ^ |g?i|, the actual 
trajectory will oscillate with time but will not intersect 
the adiabatic trajectory. As a result, the computational 
error will have some oscillations but never reach zero. 

Note that in this section we have chosen 
N = 10, M = 1 for Figs. |3|4| k|5l just for the clar¬ 
ity of the figures. The essential conclusions drawn from 
these figures are the same for larger N. 


VII. CONCLUSION 

In sum, with the aid of hierarchical theory [12], we have 
shown that it is very effective to shrink the computa¬ 
tional error by controlling the time derivative of s(t) at 
the beginning and end of the evolution. Our numerical 
results with six typical adiabatic paths show that a path 
of higher orders (smoother path by intuition) leads to er¬ 
rors of orders of magnitude smaller. Or, for an allowed 
computational error, the algorithm with a higher-order 
path can be order-of-magnitude faster. The large devia¬ 
tion from the ground state in the middle of the evolution 
is not essential as long as it does not break the adiabatic- 
ity too much. Although we have focused on quantum 
search, our results and method are general and can be 
applied to other quantum adiabatic algorithms. 

This work is supported by the NBRP of China 
(2013CB921903,2012CB921300) and the NSF of China 
(11274024,11334001,11429402). 
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Details on first order deviation 

The first order deviation from instantaneous fixed 
points ( p , q) can be written as 

p(t) = p[s(f)] + Sp, q(t) = g[s(f)] + Sq. (22) 


8 

The reason why first-order derivatives of H c do not 
appear on the right hand side of Eq.(23) is simply 
because (p, q) is the fixed point. 


Now we consider s(t) changes slowly with time and 
examine the dynamics of (5p,Sq). We have 


Eq.(23) becomes 


dp 

dt 

dq 

dt 


dp . 

d^ S 

dq . 


dSp 

dt 

dSq 

dt 


(25) 


/dSp\ 
dt 


dSq 

w 


= r„(s) 




( d /\ 

os 

dq 

\~dsJ - 


(26) 


It can be shown that the determinant | To | does not 
vanish as long as the energy levels of H s are non degen¬ 
erate. Previous work [121 HSUTBj shows that ( 5p,5q) are a 
canonical pair and Eq.(26) can be derived from following 
Hamiltonian, 


H 1 (s,s) 


1 ,d 2 H c 

2 dq 2 


~ B l ) 2 + 


( d 2 H c 

dqdp 


)q^q(Sq- B 1 )(Sp 


Ai)+ 


(27) 


First we consider s is fixed. Using Hamilton equations of 
motion and Talyor expansion to the first order of ( Sp , Sq ), 
we have 


1 ( d 2 H c 

2 dp 2 


)p,q(Sp ~ ^l ) 2 


where 


ii 

<) 

( d 2 H c 

d 2 H c \ 

dqdp 

dqdq 

d 2 H c 

d 2 H c 

V dpdp 

dpdq J 


(23) 


where (Ai,B\) is the center of first order deviation 
(Sp, Sq) and defined as 





(f\ 

OS 



(28) 


Higher order derivations of (p, q ) can be derived in similar 
way. 


(24) 





















